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Abstract. The problem of spectra formation in hydrodynamic approach to 
A+A collisions is discussed. It is analyzed in terms of the two different ob- 
jects: distribution and emission functions. We show that though the process 
of particle liberation, described by the emission function, is, usually, contin- 
uous in time, the observable spectra can be also expressed by means of the 
Landau/Cooper-Frye prescription. We argue that such an approximate duality 
results from some symmetry properties that systems in A+A collisions reach to 
the end of hydrodynamic evolution and reduction of the collision rate at post 
hydrodynamic stage. 

Keywords: hydro-kinetic approach, distribution functions, emission functions. 
PACS: 24.10.Nz, 24.10.Pa, 25.75.-q, 25.75.Gz, 25.75.Ld. 



1. Introduction 

The two main features of high energy nucleon and nuclear collisions, namely, strong 
interactions and multiparticle production, suggest an idea to Heisenberg, Watagin 
and Fermi [1] that the systems created in those collisions could be considered as 
thermal. The Landau hydrodynamic model for multi-hadron production [2] ap- 
pears as a method that is alternative to S-matrix one, where the asymptotic states 
at t = ±oo are not considered but instead the detail space-time evolution of ther- 
mal hadronic matter is described basing on the local energy-momentum conservation 
laws. Instead of the initially well defined asymptotic state of colliding particles, this 
approach uses initial conditions of hydrodynamic expansion that, in fact, strongly 
depend on pre-thermal stage of hadronic or nuclear collisions. At the RfflC ener- 
gies the pre-thermal stage is formed in collisions of the two groups of specifically 
distributed partons associated with colliding nuclei that results, apparently, in very 
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dense state of partonic matter: Color Glass Condensate (see, e.g., [3]). Current 
phenomenological and pQCD estimates give the (proper) times To of thermalization 
of that matter between 0.6 and 3 fm/c. Different phenomenological assumptions 
are used to fix the details of thermalization; well known examples are Landau initial 
condition: compressed static disk at t — t [2] and Bjorkcn one: quasi-inertial flow 
at proper time t = (t 2 - z 2 ) 1 / 2 = t [4]. 

The another important aspect of the hydrodynamic approach is the final state of 
colliding system. If one uses hydrodynamic equations till infinitely large times, the 
result will be the infinitesimally small final densities that gives no chance to restore 
microscopic (or particle) picture from continuous medium approach. According to 
the Landau criterion of freeze-out, the hydro-evolution stops and all particles be- 
come free when fluid elements reach the temperature that is equal to the mass of 
the lightest hadron (pion): T/. D . = l//3/. Q . = m^. At that temperature the mean 
free path in pion gas becomes to be equal to transverse radius of hydrodynamic 
tube in p + p(p) collisions and system decays. This qualitative criterion is in mys- 
terious qualitative agreement with hydrodynamic descriptions of very wide class of 
high energy collisions: from p + p(p) to heavy ion A+A (see, e.g., [5]) where the 
correspondent fitting freeze-out temperature turns out to be 120-^150 MeV. This 
universality seems like a puzzle since it does not take into account the real dynamics 
and system sizes. Moreover, the results of many studies of A+A collisions based 
on cascade (transport) models contradict to an idea of sudden freeze-out at some 
fixed temperature. The particles escape from the system during the whole period 
of its evolution and do not demonstrate the local equilibration at the late stages. 
Though the pure hadronic cascade models as well as hybrid "hydro + cascade" 
models [6] fail to describe properly experimental data, especially the HBT radii in 
A+A collisions, the problem of spectra formation in the hydrodynamic approach in 
itself is very serious and only stress the puzzle of the successful application of the 
hydrodynamic models to multiparticle production processes. 



2. Distribution and emission functions and observables 

To clarify the problems we will follow to the basic ideas of Ref. [7] and express the 
observables through both distribution and emission functions of expanding systems. 
We start from the Boltzmann equations (BE) as more general than the hydrody- 
namic ones. The BE for the distribution function f(x,p) in the case of no external 
forces has the form: 

^^p-=F^ in (x,p)-F loss (x,p). (1) 

The term F gam and F loss are associated with the number of particles which re- 
spectively came to the point (x,p) and leave this point because of collisions. The 
term F s (x,p) — R(x,p)f(x,p) can easily be expressed in terms of the rate of 
collisions of the particle with momentum p, R(x,p) =< av re i > n{x). The term 
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pgain jjgg more complicated integral structure and depends on the differential cross- 
section. Let us split the distribution function at each space-time point into two 
parts: f(x,p) = f lnt (x,p) + f esc (x,p), x = (i,x). The first one, f int (x,p), describes 
the fraction of the system which will continue to interact after the time t. The 
second one, f esc (x,p), describes the particles that will never interact after the time 
t. According to the probability definition 

fesc(x,p) = P(x,p)f(x,p), (2) 

where escape probability P(x,p), or probability for any given particle at x with 
momentum p not to interact any more, propagating freely, can be expressed explic- 
itly in terms of the rate of collisions along the world line of the free particle with 
momentum p through the opacity integral 

P{x,p) = exp ^- j dt'R(x',p)j (3) 

Since f esc (x,p) is formed from the particles suffering last collisions at space-time 
point x it is associated with term PF 9ain and form the emission function S [7] 

p^f eS c(x lP ) = p P(x,p)F° am (x,p) EE S(x,p). (4) 

For initially finite system with a short-range interaction among particles, the 
system becomes free, in fact, at large enough times t ou t, so P(x,p) — ► 1 and 
fesc(x,p) — > f(x,p) in this limit. Therefore, to describe the inclusive spectra of 
particles 

o dN . , . o o dN i + + \ / r \ 

P dp" = ^ ap) ' PlP2 d Pl dp 2 = ^>i>2 a Pi a P2>. ( 5 ) 

the asymptotic equality f esc (x,p) = f(x,p) can be used, replacing the total distri- 
bution function / in all irreducible averages in (5), 

(a+a P2 ) = / da f _ l p^exp{iqx)f(x 1 p), (6) 

J <?out 

by f esc . Here, p = (p 1 +p 2 )/2 , q = pi - pi and the hypersurface a out just gener- 
alizes t ou t ■ Applying the Gauss theorem and recalling that d^lp^exp^qx)] = for 
particles on mass shell, one obtains, using respectively general equations (4) and 
(1) and supposing their analytical continuation off mass shell, 

(a+a P2 ) = P» f daJ esc (x a , P y c > x + r° Ut d 4 xS(x,py> x , (7) 

(a+a P2 ) = p" f daJ(x 0l p)e^+p o rd 4 x(F^ n (x,p)-F loss (x,p))e^(8) 

Jan J&Q 
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where S(x,p) is defined through the product PF aain by Eq. (4) , f eS c{xo,p) cor- 
responds to the portion of the particles, which is already free at initial time i 0j 
or, more generally, at the initial hypersurface cto, and f(x n ,p) is the distribution 
function at er - 

Thus, the use of escaping function as the asymptotic interpolation to the solu- 
tion of BE is equivalent to taking, as the source function for the spectra and corre- 
lations, the 4- volume emission function S = p°PF gam together with direct emission 
fesc(xo,p) from an initial 3D hypersurface <jq. The Landau criterion of freeze-out 
of locally equilibrium (l.eq.) hydrodynamic momentum spectra and correspondent 
Cooper-Frye prescription, defined by Eq. (6) with substitutions a ou t — > <r/. . and 
/ —* J l.eq. j treats particle spectra as results of rapid conversion of a l.eq. hadron 
system into a gas of free particles at some hypersurface ct/. .. Formally, it corre- 
sponds to taking the cross-section tending to infinity at t < t af o (to keep system 
in l.eq. state) and zero beyond t (7f o . Then P(t,x,p) = 9(t — t afg (xj) (and so 
f esc = at t < t af o ), and S = jAPF'™ = p°S{t - t af Jx))f Leq . in Eq. (7). The 
proportionality between the S and fi. eq ., like S(x,p) = p{T)fi. eq .(x,p) is used in 
many papers devoting to a description of the data in a hydrodynamically motivated 
way (see review [8]). One should understand, however, that in the realistic case 
of no sudden freeze-out the emission function S(x,p) loses completely its propor- 
tionality to the distribution function f(x,p), it is just two different objects! This 
was shown by direct calculations in Ref. [7] and also is clear from Eqs. (2), (4) 
since the escape probability P(x,p) in finite systems is very sensitive to an asym- 
metry of the positions x as for the effective "boundary" and, thereby, the S(x,p) 
becomes anisotropic in momentum p in the rest frame of a fluid element unlike to 
the hydrodynamic distribution function fi. eq .{x,p). 

3. The duality of the hydro and kinetic approaches 

As it was advocated in the recent papers [9] the perfect hydrodynamics is a good 
approximation for the earlier stage of the matter evolution because of the big cross- 
section of the interaction among colour and white quasi-hadronic states in QGP. 
As for the matter evolution at the post hadronization stage in A+A collisions it 
was found the approximate equivalence between chemically frozen hydrodynamics of 
hadron gas and its evolution within the cascade approach in the temperature region 
above 0.12 GeV [10]. Below this region the formation of spectra is continuous in time 
with fairly long "tails" of the emission. That process is characterized by the emission 
function S(t, x, p) that is far from thermal. The similar situation was analyzed in 
Ref. [7] based on the exact analytic solution of the BE for expanding fireball. 
Unlike to very complicated structure of the emission function S(t, x, p) that was 
found in [7], the spectra, intcrfcrometry radii, averaged phase-space densities have 
simple and clear analytical forms and correspond to the Cooper-Frye prescription 
for the thermal distribution function fi. eq .(x,p) in the fireball before it starts to 
decay. The explanation is based on the duality of Eqs. (7) and (8) that express the 
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spectra in terms of either the emission or distribution functions. While the emission 
function, that is proportional to F gam , can have a significant non-zero value in wide 
space-time region, the integral of the difference F gam — F loss over this region could 
be zero, as in the example discussed in Ref. [7], or small. The latter is typical 
when the system expands in spherically symmetric way. The analysis of different 
hydrodynamic solutions demonstrates that the velocity field of expanding systems 
tends typically to a spherically symmetric one at the late stage of evolution, at least, 
in the central region where the low px spectrum forms. As for the high px spectrum 
formation, one can expect that correspondent particles are radiated mainly from 
the periphery of the system at the earlier times, because of large hydrodynamic 
velocities and fast transition to free streaming there: it was argued in Refs. [7] and 
[11]. Thereby, in a rough approximation, one can apply the generalized Cooper- 
Frye prescriptions (8) putting there <r = &f.o., taking into account the possible 
p-dependence of hypersurface <Jf. .(p) at high Vt h 0.8 — 1 GeV and neglecting the 
integral over the 4D region situated beyond of the (jf. .{p). That is some kind of 
duality in the description of the spectra and the intcrfcrometry data basing cither 
on (thermal) distribution functions fi, eq ,(x,p) which characterize the system just 
before decay begins or on the emission function S(t, x, p) that describes the process 
of continuous particle liberation during the decay of system. 

In fact, any kind of hydrodynamic or Boltzmann kinetic approaches loses its 
applicability when the particle mean free paths become compatible with lengths of 
homogeneity in the system. The universality of the Landau freeze-out temperature, 
T ~ , an independence on the concrete dynamics and correspondent homogene- 
ity lengths is due to the pions are dominated in hadron system and their mean free 
paths, l/an(T(x)), start to increase exponentially ~ exp(/3m 7r ) with j3 — 1/T(x) 
when the temperature falls down below as it follows from the analytic repre- 
sentation of the thermal density n(T). Thus the temperature T ~ m P i is just lower 
boundary of the region of applicability of hydrodynamics in wide class of nucleon 
and nuclear collisions. It does not mean that the hadrons stop to interact then at 
post hydrodynamic stage but the momentum spectra do not change significantly 
especially if the above discussed conditions (small value of F 9am , symmetry of ex- 
pansion at the late stage, etc.) are satisfied and so the integral of F 9am — F loss in 
Eq. (8) is small at that stage. 

4. Conclusions and outlook 

We show that universal Landau freeze-out temperature corresponds to lower bound- 
ary of the applicability of hydrodynamics that is similar in different collisions. The 
duality in spectra description between the (generalized) sudden freeze-out prescrip- 
tion, that utilizes the distribution functions /, and the detail picture of the particle 
liberation process, based on the emission function S, is argued. The replace of the 
complicated emission process by the simple Landau [2] criterion of sudden freeze-out 
at T « 0.12 -7- 0.15 GeV is, of course, a rather rough approximation. Nevertheless, 
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as it was discussed in Ref. [12], the momentum-energy conservation laws, peculiari- 
ties of almost isoentropic and chemically frozen evolution as well as some symmetry 
features of the late stage of hydro evolution minimize the correspondent uncertain- 
ties. Note, that the duality does not mean and even excludes the parametrization 
of emission function in the form £ ~ /;. eg . excepting for the case of real sudden 
freeze-out which is, probably, very non-realistic. 

The approximation discussed is, unfortunately, not well controlled in practical 
utilization and we hope to develop an approximate method for spectra calculations 
within the hydro-kinetic approach [7] based on the escape probabilities and gen- 
eralized relaxation time approximation. The method will combine the advantages 
of hydrodynamic approximation and microscopic (kinetic) approach. The former 
allows one to incorporate the complicated evolution of the system at the possible 
phase transitions encoded in corresponding equation of state; the latter makes pos- 
sible to evaluate the observable spectra taking into account the non-equilibrated 
character of their formation. As a result we hope coherently explain all totality of 
hadronic observables of various experiments where existing dynamic models fail to 
accommodate the majority of experimental data. 
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